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ABSTRACT 

Context. Coronal bright points (CBP) are ubiquitous small brightenings in the solar corona associated with small magnetic bipoles. 
Aims. We derive the solar differential rotation profile by tracing the motions of CBPs detected by the Atmospheric Imaging Assembly 
(AIA) instrument aboard the Solar Dynamics Observatory (SDO). We also investigate problems related to detection of coronal bright 
points resulting from instrument and detection algorithm limitations. 

Methods. To determine the positions and identification of coronal bright points we used a segmentation algorithm. A linear fit of their 
central meridian distance and latitude versus time was utilised to derive velocities. 

Results. We obtained 906 velocity measurements in a time interval of only 2 days. The differential rotation profile can be expressed 
as oJrot - (14.47 ± 0.10 + (0.6 ± 1.0) sin^(Z7) + (-4.7 ± 1.7) sin"'(Zi))°day“*. Our result is in agreement with other work and it comes 
with reasonable errors in spite of the very short time interval used. This was made possible by the higher sensitivity and resolution of 
the AIA instrument compared to similar equipment as well as high cadence. The segmentation algorithm also played a crucial role by 
detecting so many CBPs, which reduced the errors to a reasonable level. 

Conclusions. Data and methods presented in this paper show a great potential to obtain very accurate velocity profiles, both for 
rotation and meridional motion and, consequently, Reynolds stresses. The amount of coronal bright point data that could be obtained 
from this instrument should also provide a great opportunity to study changes of velocity patterns with a temporal resolution of only 
a few months. Other possibilities are studies of evolution of CBPs and proper motions of magnetic elements on the Sun. 

Key words. Sun: rotation - Sun: corona - Sun: activity 


1. Introduction 

We present a new solar rotation profile obtained by trac¬ 
ing the motions of coronal bright points (CBPs) observed by 
Atmospheric Imaging Assembly (AIA) instrume nt on board 
the S olar Dynamics Observatory (SDO) satellite (iLemen et al.l 
I2012I) . 

The most frequently used and oldes t tracers of the solar dif¬ 
ferential rotation profile are sunspots (iNewt on & Nu m il 19^ 


EIT observations in 28.4 nm channel and 


(1201 Oh used SOHO- 


Karachik et alJ (l2006h 


used 19.4 nm SOHO/EIT channel. iKarivappal (l2008h also used 
Hinode/XRT full-disk images to determine the solar rotation 
profile. 

Other tracers are used as well: mag netic fields 

l^mm et all 


(IWilcox & HowardI 119701: iSnodgrassI 
11993 ) and Ha filaments ( Braisa et alJll991l 


jnn 

]). Apar 


_ Apart from tracers, 
d ( [Howard & HarvevI 

ichlT 


Howard et al.lll984tlBalthasar et alJ[l98i^ : lBraisa et al.ll2002ah . 

One of the advantages of using sunspots is very long time cov¬ 
erage. On the other hand, there are numerous disadvantages; 
sunspots have complex and evolving structure, their distribution 
in latitude is highly non-uniform and it does not extend to higher 
solar latitudes. The number of sunspots is also highly variable 
during the solar cycle which makes measurements of solar differ¬ 
ential rotation profile almost impossible during solar minimum. 

CBPs are more uniformly distributed in latitude and are nu¬ 
merous in all phases of the solar cycle. They also extend over all 
solar latitudes. They have been used as tracers of solar rota tion 
since the beginning of the space age (iDupree & Henzell19721) . In 
recent years there are numerous studies i nvestigating solar differ¬ 
ential rotation by using CBPs as tracers. Kari^aroa ( 2008h: Hard 
(1200^ used Yohkoh/SXT d at a whilelBraisa et al.1 (l2001ll200M 
120041) : IVrsnak et alJ (l2003l) : IWohl et all 


Send offprint requests to: D. Sudar, e-mail: : davor.sudart® gmail.com 


Dopp ler measurements can also be used 
IT^ Ullrich eT^IT^ISnodgrass & Ulrichil990h . 

Helioseismic measurements also show differential rotation 
below the photosp here all the way down to the bottom o f 
the convective zone (iKosovichev et al.iri997l : [^hou et alJll998l) . 
Eurther dow n, the rotatio n profile becomes uniform for all lati¬ 
tudes (cf. eg. lHow 3120091) . 

Eor further details about solar rotation, its importance for 
solar dynamo models, and comparison of rotation measure¬ 
ments b etween d i fferent source s, see t h e reviews by Schrot^ 
(119^: [HowardI (119^: [Be^ (l2000l): lOssendriiveil (| 200 ^ 
[Rudiger & Hollerbachl (|2004 : [siiT (l2004 : iHowg (1200^ : 
iRozelot & Neineil^2009h. 

In this work we use CBP data obtained by SDO/AIA over 
only two days to assess the quality of the data, identify sources of 
errors and calculate the solar differential rotation profile. We will 
also investigate the possibility of using CBP data from SDO/AIA 
for further studies of other related phenomena (meridional flow, 
rotation velocity residuals and Reynolds stress). 

CBP data from SDO/AIA were also used in other works. 
iLorenc et al.l ( 1201 2|) discussed rotation of the solar corona based 
on 69 structures from 674 images detected in 9.4 nm channel 
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using an interactive method of detection. iDorotovic et al.l (l2014t) 
presented a hyb r id algor ithm for detection and tracking of CBPs. 
iMcIntosh et al.l (l2014b) used detection algorithm presented in 
their previous paper ( McIntosh & Gurmim l2005h to identify 
CBPs in SDO/AIA 19.4 nm channel and correlate their proper¬ 
ties with those of giant convective cells. Using more SDO/AIA 
data an d extending analysis back to SOHO era, iMcIntosh et alJ 
(12014ah concluded that CBPs almost exclusively form around 
the vertices of giant convective cells. 


2. Data and reduction methods 


We have used data from the AIA instrument on board the SDO 
satellite (iLemen et al.ll2012|) . The spatial resolution of the instru¬ 
ment is «j0.6”/pixel. For comparison, SOHO/EIT resolution is 
2.629”/pixel while Hinode/XRT has a resolution of 1.032”/pixel. 

To obtain positional information for the coronal bright points 
(CBPs), we employed a segmentation algorithm which uses the 
19.3 nm AIA channel data to search for localized, small intensity 
enhancements in the EUV compared to a smoothed background 
intensity. More details abo ut the detection algor i thm, w hich is 
similar t o the algorithm by IMcIntosh & GurmanI (l200-5li . can be 
found in iMartens et alJ (|2012|) . 

This resulted in measurements of 66842 positions of 13646 
individual CBPs covering two days (1st and 2nd of January 
2011). The time interval between two successive images was 10 
minutes. In top panel of Fig. [T]we show the distribution of de¬ 
tected CBPs and compare it to the full disk image of the Sun 
in the 19.3 nm channel obtained on 1st of January 2011 (bot¬ 
tom panel of the same Figure). In the bottom panel, white cir¬ 
cles show CBPs that were detected on one image by the seg¬ 
mentation algorithm. We can see that CBPs are scarce in active 
regions, partly because of difficulties in detecting them against 
such bright and variable backgrounds. 

The segmentation algorithm provides coordinates in pixels 
(centroids of CBPs on the image) and we converted them to he¬ 
liographic coordinates tak ing into account the cur rent solar dis¬ 
tance given in FITS hies (iRosa et al.lll995L Il998h . Positions of 
objects near the solar limb are fairly inaccurate. Limiting the 
data to +58°from the centre of the Sun or of the pro¬ 

jected solar disk remo ves this problem (cf. lStark & W6h]|ll981t 
fealthasar et aljrT986af) . 

As can be seen from Fig. |2] the calculated velocities show 
some scatter. This scatter arises because the shifts are fairly 
small at our 10 min cadence, and there can be signihcant varia¬ 
tions in brightness and structure of CBPs which inhuences cal¬ 
culation of the centroid points. Nevertheless, trends are visible, 
especially in azimuthal motion which is known to be a signif¬ 
icantly larger effect. The dominant azimuthal motion led us to 
approximate the CBP motion with a linear ht to calculate the 
velocities: 
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where ojsyn is a synodic rotational velocity, tOmer is a meridional 
angular velocity, is central meridian distance (CMD) and h, 



X [pixels] 



Fig. 1. Distribution of CBPs detected by the segmentation algo¬ 
rithm (top panel) and image of the Sun in the 19.3 nm chan¬ 
nel obtained by SDO/AIA on 1st of January 2011. White circles 
show detected CBPs on this image (bottom panel). 


is latitude of each measurement for a single CBP. We have also 
removed all CBPs which had less than 10 measurements of po¬ 
sition in order for linear hts to be more robust. This is equivalent 
to 100 minutes or about l°at the equator. To obtain the true rota¬ 
tion of CBPs on t he Sun we convert sy nodic velocities to sidereal 
using Eq. 7 from ISkokic et al.l (l2014l) . 

Trying to identify the same object on subsequent images with 
an automatic method is bound to result in some misidentifica- 
tion. The resulting velocities are usually very large and can eas¬ 
ily be removed by applying a simple velocity filter. Even the hu¬ 
man fa ctor can introduce such errors. Eor example. ISudar et^ 
(1201 4t) analysed solar rotation residuals and meridional mo¬ 
tions of sunspot groups from the Greenwich Photoheliographic 
Results and found that they had to use a filter 8< a>rot <19°day“' 
for rotational velocity in order to eliminate these erroneous mea¬ 
surements. The Greenwich Photoheliographic Results catalogue 
is being investigated and revised partly in order to remove such 
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t [days] 



Fig. 2. Motion of a single CBR Top panel: latitude, b, over time, 
f; bottom panel: central meridional distance, CMD over time, t. 


problems (IWillis et al.ll2013 JSlErwin et al.ll20T^ . In this work 
we have also used a 8< Urot <19°day"* filter for rotational ve¬ 
locities to remove such outliers. In addition, we applied a merid¬ 
ional velocity hlter of -4< Wmcr <4°day“' to remove further out¬ 
liers. 

After completing all the procedures described above, we ob¬ 
tained 90 6 velocity measurements by tra cing CBPs over just 
two davs. lOlemskov & KitchatinovI (l2005h pointed out that non- 
uniform distribution of tracers can result in false flows. This ef¬ 
fect is most notable for meridional motion and rotation veloc¬ 
ity residuals, but can easily be removed by assigning the calcu- 
lated velocity to the latitude of the hrst measuremen t of position 
(lOlemskov & KitchatinovI 120051: ISudar et ^ l2014t) . Although 
the effect is negligible for solar rotation, we nevertheless applied 
the correction in this work. 

It is important to keep in mind that even when the tracers are 
uniformly distributed over the solar surface, the distribution of 
tracers in latitude will be non-uniform. As we move from equator 
to the pole, the area of each latitude bin becomes smaller, so we 
observe progressively fewer tracers (~ cos b). 


3. Results 

In this work, we present an analysis of the motion of CBPs ob¬ 
served by the SDO/AIA instrument. For a better understanding 
of the results and the potential of future studies along these lines, 
it is very useful to analyse the accuracy and errors of the dataset. 
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Fig. 3. Errors of rotational velocity, cr{Urot), for each of the 906 
measurements. 



Fig. 4. Distribution of rotation velocity errors (ojrot) in helio¬ 
graphic coordinates. Error scale is in ° day“^ 

In Eig. [3]we show errors of the calculated rotational veloc¬ 
ities, cr{Urot), for each CBP, which resulted from errors in the 
linear htting of longitude vs time measurements. Although the 
errors can go up to 3°day“', the majority is below l°day“'. In 
Eig.|4]we show these errors in heliographic coordinates to check 
their spatial distribution on the solar surface. Larger errors are 
shown with brighter shades and we can see that these roughly 
correspond to the positions of active regions shown in the bottom 
panel of Eig. [1] This correlation with active region is probably a 
consequence of the detection algorithm design and difficulties in 
detection of CBPs over a bright, variable background. 

In Pig. |5]we show the distribution of CBPs in heliographic 
coordinates with arrows indicating the velocity vector. As ex¬ 
pected, the dominant effect is that of solar rotation. 

The latitu dinal dependence of rotational velocity is usually 
expressed as dHoward & Harve^ll970l:ISchroterlll985h : 

(^rotib) - A + B sin^ b + C sin"^ b, (3) 

where b is the latitude. Parameter A represents equatorial ve¬ 
locity, while B and C depict the deviation from rigid body ro¬ 
tation. The problem with Eq. |3] is that the functions in this 
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Fig. 5. Distribution of CBPs in heliographic coordinates with ar¬ 
rows showing the direction and strength of the velocity vector. 


Fig. 6. Solar differential rotation profile obtained with data from 
SDO/AIA. Open circles are individual measurements, while the 
solid line is the best fit defined by Eq. [3]for the A + B + C case. 



expression are not ortho gonal, so the parameters are not in¬ 
dependent of each other (Duvall & Svalgaard 1978t ^odgrassl 


11984 ISnodgrass & HowardI 119851 : ISnodgrass & Ulrichl Il99d) . 
This crosstalk among the coefficients is particularity bad for B 
and C. The effect of crosstalk does not affect the actual shape of 
the fit (urotib)), but it creates confusion when directly compar¬ 
ing coefficients from different authors or obtained by different 
indicators. 

There are various methods to alleviate this problem. 
Frequently C is set to zero since its effect is noticeable only 
at higher latitudes. This is almost a standard practise when ob¬ 
serving rotation by tracing sunspots or sunspot g roups because 
their p ositions do not extend t o high latitudes ( Howm;det^ 
1984 Balthasar et al. Il986bt Pul^inen & TuominenI 1199^ 


Brai^a et al .112002a : Sudar et ^I20l4).~ 

Another method to reduce the crosstalk problem is to set the 
CjB ratio to so me fixed value. IScherrer et al.l (Il980l) set the ratio 
CjB - I while luirich et al.l (Il988ll . after measuring the covari¬ 
ance of B and C, set the ratio to ClB = 1.0216295. 

In Fig. I^we show individual measurements of rotational ve¬ 
locities, a>rot, with respect to latitude, b, as open circles. We in¬ 
dicate, with a solid line, the best fit to the data using a func¬ 
tional form given in Eq. |3] Coefficients of the fit are given in 
Table [U We also fitted the rotation profile for Northern and 
Southern sol ar hemisphere sep arately because of possible asym¬ 
metry (cf.eg. lWohl et al.l|2010l) and show the results in the same 
table. Coefficient A shows a larger value in the Southern hemi¬ 
sphere for all 3 fit functions. I.Turdana-Sepic et aTI d201 ih re¬ 
ported that coefficient A is lar ger when solar activity is smaller . 
According to the SIDC data (ISIFSO World Data Centeiil201 Ih 
we can see that Northern hemisphere is more active both when 
looking at the monthly smoothed means and daily sunspot data, 
consistent with IJurdana-Sepic et al.ldToi ih . However, judging by 
the errors of the coefficients, the difference between South and 
North is statistically low and the hypothesis that this is a result of 
asymmetric solar activity needs to be verified with a larger data 
sample. 

In Fig. |7] we show a comparison between different fitting 
techniques: A + B + C (solid line), A + B - C (dashed line) 
and A + B,C - ^ (dotted line). In the same figure we also show 


Table 1. Coefficients of the solar rotation profile. 


Type 

A [°day *] 

B [°day *] 

C [“day-*] 

n 

A, B,C 

i4.47±0.10 

-t-0.6±1.0 

-4.7±1.7 

906 

A,B = C 

i4.59±0.07 

-1.35±0.21 

-i.35±0.21 

906 

A,B,C = 0 

i4.62±0.08 

-2.02±0.33 

0 

906 

Northern hemisphere 

A, B,C 

i4.43±0.13 

+0.8±1.5 

-5.6±3.0 

461 

A,B = C 

i4.55±0.10 

-1.35±0.35 

-i.35±0.35 

461 

A,B,C = 0 

i4.57±0.10 

-1.92±0.52 

0 

461 

Southern hemisphere 

A, B,C 

i4.50±0.15 

+0.7±1.4 

-4.8±2.3 

445 

A,B = C 

i4.65±0.1i 

-1.39±0.28 

-i.39±0.28 

445 

A,B,C = 0 

i4.69±0.12 

-2.14±0.45 

0 

445 


average values of ojrot in bins 5° wide in latitude, b, with their 
respective errors. 


4. Discussion 


In Table |2] we show a comparison of the solar differential profile 
(Eq. |3]l from a number of different sources, including the results 
from this paper (Table[T]i. Since we found no statistically signif¬ 
icant difference between Northern and Southern hemisphere we 
only include the results for both hemispheres combined. Results 
in Table|2]come from a wide variety of different techniques, trac¬ 
ers and instruments. 

ISnodgra^ (11984) suggested that the rotation profile should 
be expressed in terms of Gegenbauer polynomials since they are 
orthogonal on the disk. This eliminates the cross-talk problem 
between coefficients in Eq. |3] Using the expansion in terms of 
Gegenbauer polynomials solar rotation profile becomes: 

oJroAb) ^ AcT^isinb) + BcT^isinb) + CcTlisinb), (4) 


where Ac, Bq and Cc are coefficients of expansion and T^isinb), 
r/ l^sinfc) and r|(sinf?) are Geg enbauer polynomials as defined 

bv ISnodg rass & HowardI (1 1 98.51) in t heir equation (2). _ 

As ISnodgrass & Howardl(ll98-5h : [^odgrass & lJlrichl(ll990h 
pointed out, the relationship between coefficients A, B and C 
from standard rotation profile (Eq. |3]l and coefficients Ac, Be 
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Table 2. Comparison with some other results. 


Method/object 

time period 

A [“day 

B[°day-‘] 

C [“day-*] 

Ac [°day 

Bg [“day-‘] 

Cc [°day-‘] 

Ref 

CBPs 

1967 

14.65±0.2 



14.65 



(1) 

CBPs 

1994-1997 

14.39±0.01 

-1.91±0.10 

-2.45±0.17 

13.80 

-0.709 

-0.117 

(2) 

CBPs 

1998-1999 

14.454±0.027 

-2.22±0.07 

-2.22±0.07 

13.82 

-0.740 

-0.106 

(3) 

CBPs 

1998-2006 

14.499±0.006 

-2.54±0.06 

-0.77±0.09 

13.93 

-0.611 

-0.037 

(4) 

CBPs 

1-2 Jan 2011 

14.47±0.10 

-h0.6±1.0 

-4.7±1.7 

14.19 

-0.507 

-0.224 

(5) 

CBPs 

1-2 Jan 2011 

14.59±0.07 

-1.35±0.21 

-1.35±0.21 

14.20 

-0.450 

-0.064 

(5) 

CBPs 

1-2 Jan 2011 

14.62±0.08 

-2.02±0.33 


14.22 

-0.404 


(5) 

sunspot groups 

1853-1996 

14.531±0.003 

-2.75±0.05 


13.98 

-0.550 


(6) 

sunspot groups 

1874-1976 

14.551±0.006 

-2.87±0.06 


13.98 

-0.574 


(7) 

sunspot groups 

1878-2011 

14.499±0.005 

-2.64±0.05 


13.97 

-0.528 


(8) 

sunspot groups 

1880-1976 

14.37±0.01 

-2.59±0.16 


13.85 

-0.518 


(9) 

sunspots 

1921-1982 

14.522±0.004 

-2.84±0.04 


13.95 

-0.568 


(10) 

sunspot groups 

1921-1982 

14.393±0.010 

-2.95±0.09 


13.80 

-0.590 


(10) 

Ha filaments 

1972-1987 

14.45±0.15 

-0.11 ±0.90 

-3.69±0.90 

14.11 

-0.514 

-0.176 

(11) 

magnetic features 

1967-1980 

14.307±0.005 

-1.98±0.06 

-2.15±0.11 

13.73 

-0.683 

-0.102 

(12) 

magnetic features 

1975-1991 

14.42±0.02 

-2.00±0.13 

-2.09±0.15 

13.84 

-0.679 

-0.100 

(13) 

Doppler 

1966-1968 

13.76 

-1.74 

-2.19 

13.22 

-0.640 

-0.104 

(14) 

Doppler 

1967-1984 

14.05 

-1.49 

-2.61 

13.53 

-0.646 

-0.124 

(15) 

Helioseismology 

1996 

14.16 

-1.63 

-2.52 

13.62 

-0.662 

-0.120 

(16) 

Helioseismology 

Apr 2002 

14.04 

-1.70 

-2.49 

13.49 

-0.672 

-0.119 

(17) 


Re ferences. (1) iDunree ife Henze! (119721): 

(6) iPulkkinen & Tuo minenI 1 19981): 17) Balthar_ _ _ _ _ _ _ _ _ 

1 lilSraisa et alj Jl99lh: (12)]Snodgrassl ) ll983l) : (13) lKomm et all)ll993l) : (14) [Howard & Harvevl)ll97(]h : (15) ]Snodgrassl <1984 : (16) ISchou et all 
19981): (17)lKomm et al.l<20()3) 


__ _ (2) E ar al ([200^ : 13) iBraisa et al.l <2004: 14) IWdhletalJ <20l4: 15) this 

Balthasaretayj 1986b|); 18) ISudar et all 1 2014^]_(9)_ Bra jsa et al.l 2002a^]_n0) Howardet_^ 


paper; 

] 



Fig. 7. Comparison of three different fitting procedures for the 
solar differential rotation profile. Average values of a>rot in 5° 
bins in latitude, b, are also shown with their respective errors. 


and Cc from Eq. |4] is linear. Therefore, it is not necessary to 
recalculate the fits using Gegenbauer polynomials, we can com¬ 
pute Ac, Be and C c directly from A, B and C. We used the re¬ 
lationship given in ISnodgrass & HowardI (1 19851) (their equation 
(4)) since there seems to be a typo for a simi l ar rela tionship for 
coefficients C and Cc in ISnodgrass & Ulrichl <19901) . In Table |2] 
we also show the values of coefficients Ac, Be and Cc- 

Our rotational profile results are roughly consistent with all 
the previously published work we surveyed (Table |2l). The ac¬ 
curacy of our coefficients is lower when compared with other 
results, a consequence of our fairly small number of data points 
{n - 906). IWdhl et al.l <2010h . for example, had more than 50000 
data points spanning a time interval of 8 years. We have used 
data spanning only 2 days. It is therefore reasonable to expect 


that with AIA/SDO CBP data we could reach 50000 data points 
with only 4 months of data and achieve similar accuracy in solar 
rotation profile coefficients. This means that with AIA/SDO data 
it should be possible to measure rotation profile several times per 
year and track possible changes in solar surface differential rota¬ 
tion directly with a very simple tracer method. This is also true 
for meridional motion and Reynolds stress, both of wh ich prob¬ 
ably vary over the solar cycle (cf. e.g ISudar et al.ll20T4l) . 

5. Summary and Conclusion 

Using 19.3nm data from the SDO/AIA instrument at 10 min 
cadence we have identified a large number of CBPs, resulting 
in 906 rotation velocity measurements. We obtained a fairly 
good differential solar rotation profile in spite of the fact that 
we used data spanning only two days. The large density of data 
points in time is a result of several factors. The instrument itself 
(SDO/AIA) has better spatial resolution and is capable of high 
cadence (<5 min). For comparison, SOHO/E IT 28.4 nm channe l 
had a cadence of two images every 6 hours <Wohl et al.ll20ir)l) . 
In this work, the high cadence enabled us to track and mea¬ 
sure velocities of short lived CBPs which couldn’t be detected 
or accurately tracked by the comparatively large time interval 
between successive images in the SOHO/EIT 28.4 nm channel. 
Coupled with the fact th at short lived CBP s are more numer¬ 
ous than long lived ones <Braisa et al.ll2008l) . this resulted in a 
very high density of data points in time. High data density ne¬ 
cessitated the use of an automatic procedure to detect and track 
CBPs. The segmentation algorithm used here proved to be com¬ 
pletely adequute_Jo£_the_task 2 _assimil^_algoriffim^Ji^e else¬ 
where <McIntosh & Gurmanll2005l : iMcIntosh et alJl2014alf5) . 

The surface rotation profile and its accuracy obtained by he¬ 
lioseismology is seldom given in a form suitable for comparison 
with those obtained by tracer measurements. We can estimate 
from the number o f published s i gnifica nt digits in the results 
ISchou et"^ <1998l) : iKomm et al.l <2004 . given in Table |2l that 
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the accura cy is of the same o rder as better quality tracer mea¬ 
surements. IZaatri et ^ (l2009t) published the error for the coef- 
hcient B (see Eq.O and the value of 0.01 from that paper is in 
agreement with our estimate above. 

It is quite conceivable that errors in the differential rotation 
prohle coefficients would drop signihcantly when more data is 
used. From our analysis, we can expect to obtain 400-500 ve¬ 
locity measurements per day from CBPs using SDO/AIA. A 
time interval of 4 months seems adequate to obtain 50000 ve¬ 
locity measurements, which should be sufficient to match the 
most accurate res ults obtained by tracer methods (for example 
IWohl et all (12010b ). 

CBPs are also very good tracers since they extend to much 
higher latitudes than sunspots. They are also quite numerous in 
all phases of the solar cycle while sunspots are often absent in 
the minimum of the cycle. 

This opens up an intriguing possibility of measuring the 
solar rotation profile almost from one month to the next over 
an entire cycle. Such studies could provide new insight into 
mechanisms responsible for solar rotation. We already know 
that meridional motion exhibits some changes during the course 
of the solar cycle, and th e same is probably true for Reynolds 
stress. I^dar et al.l (l2014t) found by averaging almost 150 years 
of sunspot data that meridional motion changes slightly over the 
solar cycle and hinted that the Reynolds stresses are probably 
changing too. Here we have found a small asymmetry in rota¬ 
tion prohle for two solar hemispheres and suggested that this 
might be related to different solar activity levels in the two hemi¬ 
spheres. This needs to be verihed with a larger dataset though, 
as the difference in rotation prohles was of low statistical signif¬ 
icance. 

The planned SDO mission duration of 5-10 years will cover 
a large portion of the solar cycle which should result in enor¬ 
mous amount of velocity data to assist in the understanding of 
the nature and variation of solar rotation prohle. Having more 
detailed temporal resolution and direct results (without the need 
to average many solar cycles) could prove to be very informa¬ 
tive. 

A time interval of 10 minutes between successive images 
also offers a good opportunity to study the evolution of CBPs 
and possible effect this might have on the de tected surface veloc¬ 
ity helds. For example. IVrsnak et al.l (l2003h reported that longer 
lasting CBPs show different results than short-lived CBPs. 

Based on the promising results here, we will use larger 
datasets to further exploit the potential of SDO/AIA CBP data 
to determine meridional motions, rotation velocity residuals, 
Reynolds stresses and proper motions in subsequent papers. 
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